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摘要 无 碰撞 磁场 重 联 作为 一 种 将 磁 能 有 效 转 化 为 等 离子 体 动能 和 热能 的 机 制 , 已 经 被 广泛 应 用 于 解释 太阳 
炮 斑 、 地 球 磁 暴 等 各 类 等 离子 体 的 爆发 活动 . 然而 , 在 无 碰撞 重 联 区 中 反常 电阻 的 微观 物理 机 制 仍然 是 尚未 解 
决 的 基本 问题 . 在 众多 反常 电阻 的 形成 机 制 中 , 基于 磁 零 点 附近 粒子 轨道 混沌 性 产生 的 混沌 感 
最 普遍 流行 的 形成 机 制 , 但 它 的 微观 物理 图 像 却 是 最 为 清晰 的 . 回顾 了 无 碰撞 重 联 区 中 混沌 感应 电阻 的 早期 有 
究 和 基本 理论 模型 , 介绍 了 关于 混沌 感应 电阻 研究 的 新 进展 并 阐述 了 混沌 感应 电阻 未 来 的 研究 方向 . 
关键 词 磁场 重 联 , 反常 电阻 , 混沌 感应 电阻 , 太阳 : 日 晃 物 质 抛射 
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1 引言 故 磁场 重 联 几乎 不 会 发 生 . 那么 在 无 碰撞 磁场 重 联 
磁场 重 联 的 锥 形 最 初 由 Giovanelli 在 1946 年 提 © 宁 是 什么 打破 了 磁 冻 结 条 件 使 得 磁场 重 联 发 生 的 ? 

EO, 用 来 解释 太阳 耀 斑 的 产生 机 制 . 目前 , 磁场 基于 此 , 研究 者 们 提出 了 反常 电阻 6 的 概念 . 

重 联 冲 作为 解释 太阳 光斑 名、 日 园 物 质 抛射 外 、 磁 a 

暴 回 、 磁 层 亚 暴 回 以 及 极光 四 过 程 中 磁场 能 量 快速 BOOM FATALE, 研究 者 发 现 该 模型 所 

释放 的 有 效 机 制 之 一 , 已 经 被 广泛 地 研究 和 应 用 . ”项 测 的 重 联 率 很 小 , 重 联 过 程 所 需要 的 时 间 很 长, 


由 磁场 重 联 的 概念 &-9, 可 知 磁场 重 联 发 生 的 主要 ”例如 在 日 园 抛 射 喷 流 过 程 中 , 理论 上 预测 的 重 联 时 
= 特点 是 磁场 拓扑 结构 发 生 改 变 和 等 离子 体 被 加 速 间 是 月 的 量 级 , 而 典型 观测 中 的 日 冕 喷 流 时 间 却 是 
“= 和 加 热 , 其 中 磁场 拓扑 结构 发 生 改 变 一 定 是 打破 了 ” 分 钟 的 量 级 . 因此 , Sweet-Parker 模 型 常 被 称 为 | 
WER EER EO, 在 经 典 等 离子 体 理 论 中 , 等 离子 体 联 . 之 后 关于 磁场 重 联 的 工作 大 部 分 围绕 着 如 何 提 
的 电阻 是 由 于 带电 粒子 之 间 存 在 库仑 人 碰撞 产生 的 ， 重 联 率 来 开展 的 , 其 中 一 些 工 作 将 Sweet-Parker 
称 之 为 经 典 斯 皮 策 电阻 或 经 典 碰撞 电阻 11, 此 电阻 E 联 模型 中 的 经 典 磁 撞 电阻 用 反常 电阻 代替， 以 
会 破坏 磁 冻 结 条件 , 使 得 磁场 重 联 发 生 . 但 是 在 空 此 来 提高 重 联 率 05 17) 此 外 , Petschek 在 1964 提 出 
间 等 离子 体 中 , 等 离子 体 参 量变 化 的 特征 时 间 要 远 了 Petschek 重 联 模型 3 , 在 该 模型 中 磁 能 不 仅 能 通 
远 小 于 带电 粒子 之 间 的 平均 碰撞 时 间 , 因此 由 库仑 过 欧姆 加 热 来 耗 散 , 还 能 通过 慢 激 波 将 磁 能 转化 为 
磁 撞 产生 的 电阻 非常 小 且 不 足以 破坏 磁 冻 结 条 件 ， 等 离子 体 的 动能 和 热能 , 此 模型 被 称 为 快 重 联 . 但 
2022-09-23 收 到 原稿 , 2022-12-21 收 到 修改 稿 
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Fe, 数值 模拟 研究 发 现 Petschek 重 联 模型 是 不 稳定 S 随 时间 t 的 变化 , 右 图 是 混沌 感应 电阻 ( 实 线 ) 与 经 
的 89. 随后 的 研究 发 现 快 重 联 稳定 发 生 需要 不 均匀 典 碰撞 电阻 (虚线 ) 随 平均 自由 程 的 变化 . Yoshida 
EPERE, 综 上 , 两 种 经 典 磁场 重 联 模型 中 等 Bd 在 1998 年 指出 磁 零 点 附近 磁场 的 强烈 变化 会 
都 需要 除 经 典 碰 撞 电 阻 以 外 反常 电阻 的 存在 , FFB 破坏 磁 矩 守恒 , 引起 了 带电 粒子 的 混沌 运动 , 从 而 
该 反常 电阻 在 重 联 区 中 分 布 应 是 不 均匀 的 . 引入 了 新 的 自由 度 , 最 终 导致 粹 的 增加 , 如 图 1 中 无 
无 碰撞 磁场 重 联 中 反常 电阻 的 本 质 是 重 联 电 侧 子 图 所 示 . 此 外 , 他 们 计算 了 不 同 平均 自由 程 下 
流 片 中 载 流 带电 粒子 运动 方向 的 随机 化 过 程 , 此 过 由 带电 粒子 混沌 运动 产生 的 混沌 感应 电阻 和 由 库 
程 使 宏观 有 序 的 能 量 转化 为 微观 无 序 的 能 量 . 使 仑 碰撞 产生 的 经 典 碰 撞 电 阻 , 发 现 有 效 电阻 要 比 经 
] 清 动 等 离子 体 波 的 波 粒 相互 作用 来 研究 反常 电 典 碰撞 电阻 高 1-2 个 量 级 , 如 图 1 中 右 侧 子 图 所 示 . 
阻 产 生机 制 是 最 普遍 的 , 例如 离子 声 滑动 波 B 忆 -2g、 这 是 首次 提出 可 以 用 无 磁 撞 磁场 重 联 中 粒子 轨道 


动力 学 阿尔 文 波 P7-30、 低 混杂 漂移 波 B-33 和 哨 声 混沌 性 产生 的 混沌 感应 电阻 解释 反常 电阻 的 产生 . 

BA. 真实 的 磁 能 耗 散发 生 在 等 离子 体 的 动力 学 Numata 等 人 发 现 只 有 在 开放 系统 中 能 量 耗 散 

尺度 上 , 因此 单个 带电 粒子 的 运动 在 磁 能 耗 散 中 才 会 持续 发 生 , 因此 他 们 研究 了 开放 系统 下 无 碰撞 
起 着 重要 作用 . 基于 磁 零 点 附近 强 磁 场 梯度 的 非 磁场 重 联 区 中 的 混沌 感应 电阻 , 并 使 用 试验 粒子 
> 线性 会 使 带电 粒子 运动 轨迹 变 为 混沌 的 B59, Yoshi- 模拟 和 系 综 统 计 的 方法 定量 计算 了 混沌 感应 电阻 


N da 等 B9 在 1998 年 提出 一 种 新 的 无 碰撞 反常 电阻 的 的 大 小 B949. 在 数值 模拟 中 , Numata 等 人 使 用 李 
z= 产生 机 制 , 指出 带电 粒子 的 混沌 运动 会 导致 载 流 粒 ” 雅 普 诺 夫 指 数 和 局 地 最 大 李 雅 普 诺 夫 指 数 分 别 衡 
© 子 定 向 运动 的 随机 化 . 该 机 制 的 核心 是 带电 粒子 在 量 单个 粒子 和 多 个 粒子 运动 的 混沌 程度 , 图 2 的 左 
O 。 不 均匀 电磁 场 中 的 混沌 运动 会 产生 等 效 的 “页 撞 *"， ” 侧 子 图 显示 了 单个 粒子 的 李 雅 普 诺 夫 指 数 (x(t) = 
~N 从 而 产生 了 反常 电阻 . 这 种 由 粒子 轨道 混沌 性 产生 (1/At) In(|d(t + At)|/|6a(t)|), daze AA ER PH Hr FEE 
的 电阻 被 称 为 混沌 感应 电阻 子 轨迹 的 距离 ) 随 时 间 的 变化 , 右 侧 子 图 显示 了 不 
本 文 主要 介绍 用 混沌 感应 电阻 解释 无 碰 接 重 同 混沌 区 半径 (图 中 局 中 多 个 粒子 的 局 地 最 大 李 
联 区 中 反常 电阻 产生 机 制 的 早期 研究 和 近期 进展 . 雅 普 诺 夫 指 数 (X(t) =< x(t) >) 随 时 间 的 变化 . 选 
第 2 节 回 顾 混沌 感应 电阻 的 早期 研究 ; 第 3 节 给 出 了 取 变 化 稳定 的 局 地 最 大 李 雅 普 诺 夫 指 数 所 对 应 的 
混沌 感应 电阻 的 基本 理论 模型 ; 第 4 节 介 绍 了 混沌 。 混沌 区 半径 , 统计 此 半径 内 带电 粒子 在 加 速 电场 
感应 电阻 在 有 限 宽度 电流 片 和 岛 链 电 流 片 中 的 最 方向 上 的 平均 速度 过 (t) 和 混沌 区 内 带电 粒子 的 个 
新 研究 进展 B7 OS), 第 5 节 给 出 总 结 并 提出 混沌 感应 Ant), 并 使 用 一 元 一 次 函数 ( 坊 (t) = at’, a 是 平 
E 阻 的 未 来 研究 方向 . 均 加 速度 ) 和 单 指数 函数 (n(t) = noexp( 一 Bt), no 
T e 是 初始 带电 粒子 个 数 , 6 是 相对 逃逸 率 ) 拟 合 , 从 而 
了 2 混沌 感应 电阻 的 早期 研究 ere a 
无 碰撞 磁场 重 联 中 反常 电阻 的 本 质 是 磁场 重 度 和 带电 粒子 离开 混沌 区 域 的 相对 逃逸 率 . 理论 上 
联 区 电流 片 中 载 流 带 电 粒 子 的 运动 方向 产生 了 在 系统 达到 平衡 状态 时 , 求 出 加 速 电场 方向 的 饱和 
随机 化 , 从 而 使 宏观 有 序 能 量变 为 微观 无 序 能 量 ， 速度 , 再 结合 由 耗 散 方程 得 到 的 饱和 速度 , 即 可 以 
故 单个 粒子 的 运动 对 反常 电阻 的 产生 至 关 重 要 . 给 出 有 效 碰撞 频率 , 从 而 计算 出 特定 等 离子 体 环境 
Schmidt 在 1962 年 发 现在 磁场 空间 分 布 不 均匀 的 磁 下 的 混沌 感应 电阻 . 他 们 使 用 日 晃 中 等 离子 体 的 典 
零点 附近 , 由 于 强人 磁场 梯度 非 线 性 的 存在 , 使 得 带 型 参数 分 别 计算 了 X 型 磁场 位 形 下 的 混沌 感应 电 
粒子 的 运动 轨迹 变 为 无 序 、 混 沌 的 B5. 图 1 的 左 和 经 典 碰撞 电阻, 发 现 混沌 感应 电阻 比 经 典 碰撞 电 
图 是 混沌 运动 ( 实 线 ) 和 周期 运动 (虚线 ) 中 动力 学 烂 。” ” 阻 高 4 个 量 级 . 


一 


m 


一 


37-2 


ChinaXiv 合 作 期 刊 


64 卷 E 桢 等 : 无 碰撞 重 联 区 中 的 混沌 感应 电阻 3H 


& T 
a 


resistivity (ohm-m) 


0 0 200 400 600 800 1000 0.0001 2 4 6 8 10 12 14 16 18 
t mean free path (m) 
图 1 Ye pelt zs oy A E hy gH ST GRAE (A ) LA YR i R EL Sg 28 TG EL BSA 程 的 变化 ( 右 )[a9] 


Fig.1 Change of kinetic entropy versus time in chaotic and periodic motion (left) and variation of chaos-induced resistivity and 


classical collision resistivity versus the mean free path (right) ° 
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四 2 单个 粒子 的 李 雅 普 诺 夫 指数 ( 左 )[89] 以 及 不 同 混沌 区 半径 中 多 个 粒子 的 局 地 最 大 李 雅 普 诺 夫 指数 随时 间 的 变化 ( 右 )I9l 


Fig.2 The Lyapunov exponent of the single particle (left)[a9] and the change of the local maximum Lyapunov exponent of 


multiple particles in different chaotic region radii (right) versus time[40 


= 地 球 磁 尾 处 的 磁场 拓扑 结构 经 常 为 X 线 电流 片 级 . 
= MOR Ha 4), 在 2014 年 , Andriyas 等 人 用 粒 磁场 重 联 可 以 分 为 完全 磁场 重 联 和 部 分 磁场 
子 轨道 混沌 性 研究 了 磁 尾 中 质子 和 氧 离 子 在 这 两 重 联 , 前 者 是 指 两 个 严格 反 平行 磁场 的 重 联 , 重 联 
种 磁场 位 形 下 的 混沌 感应 电阻 (9. 文中 使 用 了 双 指 点 的 磁场 完全 抵消 , 是 一 个 严格 的 磁 零 点 . 不 过 , 在 
BLA B(n(t!) = E; nj exp(—Byt’), j = 1,2) 拟 合 带 空间 等 离子 体 中 实际 发 生 的 磁场 重 联 大 都 是 部 分 
BE 粒子 的 相对 逃逸 率 , 并 分 析 了 庞 加 莱 截 面 和 最 大 人 磁场 重 联 , 即 其 中 重 联 点 的 磁场 并 不 是 严格 反 平 行 

李 雅 普 诺 夫 指 数 , 结果 显示 分 岔 电流 片 和 X 线 电流 的 , 因而 不 会 完全 抵消 , 重 联 后 还 有 一 个 剩余 的 分 
片 中 的 最 大 李 雅 普 诺 夫 指数 分 别 为 0.34 和 0.25, 这 量 磁场 , RA S| Syl). 在 2017 年 , Shang!" 
说 明 分 岔 电流 片 中 带电 粒子 的 运动 比 X 线 中 的 更 加 研究 了 X 型 磁场 位 形 下 引导 场 对 混沌 感应 电阻 的 
混沌 . 基于 磁 尾 等 离子 体 的 典型 参数 , 他 们 计算 了 影响 . 图 3 的 第 1 行 显示 了 3 个 带电 粒子 的 运动 轨道 
分 岔 电流 片 中 的 混沌 感应 电阻 和 经 典 磁 撞 电 阻 , 发 (Ro 为 磁场 变化 的 特征 长 度 ), 即 振荡 粒子 轨道 、 磁 
现 混沌 感应 电阻 比 经 典 碰 撞 电 阻 大 了 9-10 个 数量 漂移 粒子 轨道 和 混沌 粒子 轨道 , 第 2 行 显 示 了 这 3 个 
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粒子 所 对 应 的 磁 矩 ( 粗 线 , jv) 和 粒子 与 磁 零 点 间 的 此 外 , Shang 等 史 还 研究 了 多 粒子 系统 下 混沌 
距离 ( 细 线 ) 随 时 间 的 变化 (zs 是 特征 时 间 ), 第 3 行 显 感应 电阻 随 引 导 场 强度 的 变化 , 图 4 显示 了 不 同 加 


dl 


化 . 从 图 


中 发 现 


示 了 这 3 个 粒子 对 应 的 李 雅 普 诺 夫 指数 随时 间 的 变速 电场 下 平均 加 速度 、 相 对 逃逸 率 和 有 效 碰撞 频 
恨 荡 粒 子 轨道 对 应 的 李 雅 普 诺 夫 指 。” 率 随 引 导 场 的 变化 . 从 图 中 发 现 当 引 导 场 是 背景 磁 


数 一 直 衰减 到 0.05 以 下 , 磁 漂 移 粒 子 轨道 对 应 的 李 。” 场 的 一 半 (B, = 0.5Bo) 时 , 无 碰撞 重 联 区 中 的 有 效 
普 诺 夫 指数 稳定 在 0.05 附 近 , 而 混沌 粒子 轨道 对 ”碰撞 频率 最 大 , 即 此 时 对 应 的 混沌 感应 电阻 也 最 


应 的 李 雅 普 诺 夫 指数 稳定 在 0.1 以 上 . 基于 此 , 他 们 大 


得 到 了 混沌 粒子 轨道 对 混沌 感应 电阻 的 贡献 最 大 . 
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(i) : mirror-oscillation orbits (ii) :magnetized-drift orbits (iii) : chaotic orbits 


图 3 第 1 行 是 3 种 典 


型 运动 轨道 ; 第 2 行 是 每 


种 轨道 所 对 应 的 粒子 磁 答 ( 粗 线 ) 和 粒子 与 磁 零 点 间距 离 ( 细 线 ) 随 时 间 的 变化 ; 第 3 行 是 每 种 轨道 所 对 应 粒 
子 的 李 雅 普 诺 夫 指 数 随时 间 的 变化 [9]. 


Fig.3 The first row: three typical orbits; the second row: the magnetic moment of the particle (thick line) and the distance 


between the particle and the magnetic zero (thin line) versus time for each orbit; the third row: the Lyapunov exponent of the 


particle corresponding to each orbit versus time!*®) . 
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Ra ”不同 加 速 电 场 下 , 平均 加 速度 ( 左 )、 相 对 逃逸 率 (中 ) 以 及 有 效 碰撞 频率 ( 右 , Vere) BES] EIA HELMS 


Fig.4 The average acceleration (left panel), the relative escape rate (middle panel), and the effective collision frequency (right 


panel, Verf) versus the guide field in the different acceleration electric field[45] 


3 混沌 感应 电阻 的 基本 理论 模型 

为 了 定量 计算 带电 粒子 混沌 运动 产生 的 混沌 
感应 电阻 的 大 小 , 需 给 定 多 个 带电 粒子 , 通过 试验 
粒子 模拟 和 系 综 统 计 来 研究 混沌 感应 电阻 . 在 模拟 
中 , 选择 磁场 B = Bo(bz, by, 9) 其 中 bs = B,/Bo, 
by = B,/Bo, 5 是 一 个 常量 代表 恒定 引导 场 的 相对 
强度 , x-y 平 面 为 磁场 重 联 平面 , z 为 引导 场 的 方向 ， 
1B,、B,、B, 分 别 表示 重 联 平面 和 引导 场 的 磁场 
强度 . 沿 引 导 场 方向 上 给 出 均匀 恒定 的 加 速 电 场 ， 
即 五 = (0, 0, 五 :). 在 模拟 中 忽略 带电 粒子 之 间 的 碰 


SY 


撞 和 和 集体 相互 作用 , 则 单个 带电 粒子 的 运动 方程 为 
(1) 
u mE +v x B), 


其 中 , > 和 v 分 别 代 表 带 电 粒 子 的 位 置 和 速度 , qlm 
分 别 是 带 带电 粒子 的 电荷 量 和 质量 . 为 了 数值 模拟 的 
简便 , 使 用 下 列 的 归 一 化 参数 去 归 一 化 运动 方程 、 
磁场 和 电场 : 


B' = B/Bo, r =r/Ro, E =t TAs 

Ta = Ro/va, v = v/v , E, = E, /va Bo = Ma, 
(2) 

其 中 , Bo 为 重 联 人 磁场 的 特征 强度 , vA 是 Alfvén 速 度 ， 


MA 是 加 速 电场 的 相对 强度 . 在 模拟 中 采用 Ro = 
Va/We, 这 里 的 we 是 电磁 场 中 带电 粒子 的 回旋 频率 . 
使 用 上 述 归 一 化 参数 , 得 到 归 一 化 后 带电 粒子 


的 运动 方程 为 
dr’ i 
dt! 
a RS k 
Ww 
HEE AA BSP IAB’ = (br, by, O)AIE’ 
= (0,0, Ma). 
在 试验 粒子 模拟 中 , 可 得 到 混沌 区 域内 带电 粒 


子 在 z 方 向 的 平均 运 去 动 速度 和 带电 粒子 在 混沌 区 内 
的 剩余 个 数 , 分 别 使 用 一 元 一 次 函数 和 单 指数 函数 
来 拟 合 混沌 区 内 粒子 在 z 方 向 的 平均 运动 速度 和 带 
电 粒 子 个 数 随 时 间 的 变化 . 图 5 显示 了 带电 粒子 分 
别 在 电场 方向 上 平均 速度 和 混沌 区 域内 带电 粒子 
个 数 随时 间 的 变化 , 其 中 上 标 表 示 的 是 归 一 化 , 是 
参考 文献 40] 中 的 写法 , 与 加 ' 表 示 归 一 化 意义 一 相 


tk 


V(t) = at’, (4) 
和 

n(t') = no exp (一 Bt ) , (5) 
得 到 混沌 区 域 带电 粒子 沿 加 速 电 场 方向 的 加 速度 
a 和 带电 粒子 离开 混沌 区 域 的 逃逸 率 B. 其 中 , L(t’) 
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和 mn(t') 分 别 是 混沌 区 域内 带电 粒子 在 z 方 向 的 平均 
速度 和 区 域内 带电 粒子 的 数 密度 . 

在 理论 上 , 通过 不 断 补 充 新 的 速度 为 零 的 带电 
粒子 到 混沌 区 使 磁场 重 联 达到 稳 态 平衡 . 补充 新 的 
速度 为 零 的 粒子 之 后 , 混沌 区 内 带电 粒子 总 的 平均 


速度 为 
及 的 = 
k v,(t)noB exp(—Gt)dt + no exp(— Bt )U,(t') 
no i 
(6) 
化 简 可 以 得 到 
= PE) = F- exp (— Be). (7) 
2 HE > co 时 , 得 到 系统 达到 稳 态 时 的 饱和 速度 
- A 
a a : 
à V! = V! (tleo = =. 8 
s = V(t) 3 (8) 
0.4 
0.35 
0.3 
n 0.25 
Ja 
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同时 可 以 根据 耗 散 方程 
dV (t) 
dt 


来 计算 饱和 速度 , 其 中 wa = ver /wc 是 无 量 纲 的 有 
效 碰撞 频率 . 当 d/dt 一 0 时 , 可 得 到 饱和 速度 为 


= MA — vV; (t) (9) 


Ma 
Vs = Vi@)lajatso = 7> 
eff 


(10) 


根据 (8) 式 和 (10) 式 中 的 两 个 饱和 速度 , 可 得 到 
有 效 碰撞 频率 为 


Ma. (11) 


1/ 
Vef 一 


B 
Q 


因此 , 特定 等 离子 体 环境 下 的 混沌 感应 电阻 可 表示 
为 


T Bo 1 
eff = — Ver = — ng- 12 
Neff noe? ff ff ( ) 


Numerical Result — 


n(D=nmexp(-BD 


Number of Particles 


00 
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图 5 带电 粒子 在 电场 方向 上 平均 速度 ( 左 ) 以 及 混沌 区 域内 带电 粒子 个 数 ( 右 ) 随 时 间 的 变化 [ol 


Fig.5 The average velocity of charged particles in the direction of the electric field (left) and the number of charged particles in 


the chaotic region (right 


4 混沌 感应 电阻 研究 的 新 进展 
4.1 ”有 限 宽度 电流 片 中 的 混沌 感应 电阻 

在 研究 有 限 宽度 电流 片 中 的 混沌 感应 电阻 时 ， 
归 一 化 后 的 磁场 位 形 B' 表 达 式 为 


B= 


[40] 


versus time 


f [tanh (y!), tanh (2’ —1/2), 6], 2" > 1/2, 
< 


[tanh (y’), tanh (z + 1/2), 5], 2’ < —1/2, 
\ [tanh (y’) ,0, ô], lz'| < 1/2, 
(13) 


EP, z-y 平 面 为 重 联 平面 , > 方向 为 引导 场 的 方向 ， 
l 是 两 Y 型 中 性 点 之 间 的 距离 . 当 ! = 0 时 , 重 联 场 
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为 X 型 磁场 位 形 ; 当 ! > co 时 , 重 联 场 为 一 维 的 霍 尔 “， 率 8 基本 保持 不 变 , 故 有 效 碰撞 频率 zr 随 混沌 区 半 
BR; 当 0 < ! < co 时 , 重 联 场 中 存在 一 个 有 限 宽 = ER. Ro MNT WR); 而 磁化 区 中 平均 加 速度 a 和 
度 为 1Ro 的 电流 片 , 此 时 的 重 联 磁场 被 称 为 双 Y 型 磁 。 ”相对 逃逸 率 6 都 随 着 混沌 区 半径 有 ./ Ro 的 增加 而 减 
场 位 形 . 小 , 故 有 效 碰撞 频率 zs 随 着 混沌 区 半径 Re/Ro 增 
在 数值 模拟 中 , 使 用 2 x 10 个 带电 粒子 . Sl 一 加 而 增 大 . 此 外 , 非 磁化 区 R。 < Ro 和 磁化 区 R。 > 
0 时 , 对 应 了 X 型 磁场 位 形 , 带电 粒子 的 初始 位 置 均 。 “及 中 的 混沌 感应 电阻 都 要 高 于 转化 区 及 ~ Roth 
匀 分 布 在 一 1.0 < (zx',y) < 1.0, z = 0 中 , 且 带 电 的 混沌 感应 电阻 
粒子 具有 各 向 同性 的 麦克 斯 韦 速度 分 布 EI f(v) = TT E 
(1/(VTvur))? exp(—v?/v2), 其 中 vr = 0.3vA. 图 6 ae 
示 了 X 型 磁场 位 形 下 X 型 混沌 区 域内 平均 加 速度 粒子 的 初始 位 置 均匀 分 布 在 -5.0 < ios < ie 
on HAER GAA I igh EMME -10< < LO 2! = OF, 且 带 电 粒 了 的 初始 
区 半径 Re/Ro 的 变化 . 根据 混沌 区 域 的 大 小 , 将 其 。 EEA A TP RTE ee ISR A. 图 7 显示 了 
分 为 非 磁 化 区 R。 < Ro、 转 化 区 R。 ~ Ro 和 磁化 区 双 Y 型 磁场 位 形 和 3 种 不 同类 型 的 混沌 区 域 , 如 图 中 
— R> Ro. 从 图 6 中 , 发 现 非 磁化 区 中 平均 加 速度 w ”阴影 区 域 所 示 , 从 上 到 下 的 3 种 混沌 区 域 分 别称 之 
全 ”会 随 混沌 区 半径 R/Ro 的 增 大 而 增加 , 而 相对 逃逸 ”为 Y 型 混沌 区 域 、S 型 混沌 区 域 和 O 型 混沌 区 域 . 
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Fig.6 The average acceleration a, the relative escape rate 8, and the effective collision frequency tf versus the chaos region 


radius R./ Ro for the fixed guide field B, = 0.5Bo in the X-type magnetic configuration?) 


= 图 8 显示 了 双 Y 型 磁场 位 形 下 的 3 种 混沌 区 中 半径 R./ Ro 的 增加 而 增 大 , 而 在 O 型 混沌 区 域 却 随 

T 平均 加 速度 w。 相 对 逃逸 率 和 有 效 碰撞 频率 wr 随 Re/ Ro 的 增加 而 减 小 , 这 说 明 在 重 联 电 流 片 中 不 同 

= 混沌 区 域 半径 R./ Ro 的 变化 , 其 中 上 中 下 3 行 分 别 类 型 混沌 区 中 的 混沌 感应 电阻 存在 着 明显 差异 . 
对 应 了 图 7 中 Y 型 、S 型 和 0O 型 混沌 区 域 的 结果 . 从 选取 典型 的 光斑 前 日 园 等 离子 体 参 数 , 磁场 
图 中 发 现 平 均 加 速度 a 和 相对 逃逸 率 6 主 要 是 由 加 Bo ~ 100 Gs、 等 离子 体 密 度 no ~ 1016 m-3 和 温度 
速 电场 MA 决定 , 这 与 Shang 等 嗓 所 得 结果 一 致 , 且 T ~ 500 eV. #1 AKA, “YAY, “SAHO 


平均 加 速度 a 和 相对 逃逸 率 6 在 磁化 区 随 混沌 半 ”型 "来 分 别 表 示 X 型 泥 沌 区 域 、Y 型 混沌 区 域 、S 弄 
径 Re/ Ro 的 增 大 而 减 小 , 而 在 非 磁化 区 受 混沌 区 半 ”混沌 区 域 和 0 型 混沌 区 域 , 并 给 出 了 非 磁 化 和 碰 
径 R。/ Ro 的 影响 较 小 , 主要 的 原因 是 重 联 场 影响 了 ” 化 情况 下 有 效 碰撞 频率 的 权重 平均 值 , 根据 (12) 式 
带电 粒子 的 回旋 运动 , 从 而 阻碍 带电 粒子 的 加 速 和 ”可 计算 出 混沌 感应 电阻 mor 的 大 小 . 经 典 碰撞 电阻 
逃逸 . 此 外 , 图 中 还 显示 了 有 效 碰撞 频率 vg 在 3 种 。 ”由 公式 moa = 5.2 x 10-52 InA/T32[eV] (9 .m) 给 
混沌 区 域内 随 混沌 区 域 半径 R./ Ro 的 变化 是 不 均 。 ”出 , 其 中 2 是 离子 电荷 数 , 库仑 参数 取 In A ~ 10. 从 
匀 的 , 有 效 碰撞 频率 vg 在 Y 型 混沌 区 域 随 混 沌 区 表 1 中 发 现在 耀 班 前 日 冕 等 离子 体 中 产生 的 混沌 感 


一 


jami 


64 卷 天 文 学 报 


应 电阻 要 比 经 典 碰撞 电阻 高 5-7 个 数量 级 , 说 明 混 
沌 感应 电阻 在 无 碰撞 重 联 区 确实 可 以 产生 有 效 的 
耗 散 . 并 且 在 X 型 混沌 区 域 的 混沌 感应 电阻 要 高 于 
在 O 型 混沌 区 域内 的 混沌 感应 电阻 1-2 个 量 级 ,说 
明 混 沌 感应 电阻 在 空间 分 布 上 存在 不 均匀 性 . 


4.2 ” 岛 链 电 流 片 中 的 混沌 感应 电阻 

许多 理论 分 析 和 数值 模拟 揭示 有 限 宽度 重 联 
E 流 片 在 电阻 撕 裂 模 不 稳定 性 的 作用 下 经 常 是 不 
稳定 的 , 可 演化 为 X 型 中 性 点 和 O 型 中 性 点 交替 出 
现 的 岛 链 结构 B 5. 研究 岛 链 电流 片 中 X 型 中 性 点 
和 O 型 中 性 点 附近 的 混沌 感应 电阻 , 给 定 归 一 化 的 
磁场 位 形 芭 55 为 


= 


一 


以 


sinh(y’ /2) 
Ecos(Z//2) + cosh(y’/2)’ 

esin(x’/2) 
scos(Z'/2) + cosh(y’ /2) 
其 中 , 长 度 用 粒子 惯性 长 RA = vAy/we 来 归 一 化 , 磁 
岛 参 数 e 决 定 了 岛 链 结构 中 磁 岛 的 强度 和 位 形 . 在 
这 里 , 粒子 惯性 长 RA 可 表示 磁场 不 均匀 的 尺度 . 取 
e = 0.4 得 到 如 图 9 所 示 的 磁场 位 形 , 其 中 上 下 子 图 
中 的 阴影 区 分 别 表示 了 包含 X 型 中 性 点 的 X 型 区 域 

和 包含 O 型 中 性 点 的 O 型 区 域 . 


B= 


,0| ， (14) 


图 7 


Fig. 


3 H 
<>... nd 
Saa a 
VV Ra f 
J 下 SA 
-6 -4 -2 0 2 4 6 

a EE 
Ne oe E T 
NI A) 
—6 一 4 一 2 0 2 4 6 
a i rT A 
NMA OO oe eee AAA 
S\N TAI 
NNN 7 7 7 eee Gah 
Vo. > ae — " Na an 
-6 一 4 一 2 0 2 4 6 
x! 
双 Y 型 磁场 位 形 下 的 3 种 混沌 区 域 , 从 上 到 下 分 别 为 Y 型 混沌 区 


域 、S 型 混沌 区 域 和 O 型 混沌 区 域 371. 


7 The three chaos regions under the double Y-type 


magnetic configuration are Y-type, S-type and O-type chaos 


regions from top to bottom, respectively!” 


表 1 加速 电场 分 别 为 WMA = 0.0001、0.0005 和 0.001 时 平均 有 效 碰 撞 频 率 双 和 混沌 感应 电阻 与 经 典 碰撞 电阻 的 比 


Bnet /Neal 


[37] 


Table 1 For given acceleration electric field Ma = 0.0001, 0.0005 and 0.001, the average effective collision 


frequency V/s and the ratio of the chaos-induced resistivity to the classical Spitzer resistivity mes /ncal 
y Vef y 


are shown in the table!” 
ii Dost Nett /Meal 
0.0001 0.0005 0.001 0.0001 0.0005 0.001 
unmagnetized 0.01789 0.02585 0.03835 2.41 x 10° 3.47x10° 5.15 x 10° 
m magnetized 0.04105 0.05504 0.11460 5.51 x106 7.39x10° 1.54 x 10" 
unmagnetized 0.00205 0.00624 0.01175 2.75x10° 8.38 x105 1.58 x 10° 
etd magnetized 0.00234 0.00637 0.01594 3.14x10° 8.55x105 2.14 x 10° 
unmagnetized 0.00146 0.00427 0.00906 1.96x10° 5.73x10° 1.22 x 10° 
vas magnetized 0.00206 0.00461 0.00999 2.77x10° 6.19x105 1.34 x 10° 
unmagnetized 0.00110 0.00324 0.00784 1.47x10° 4.35x10° 1.05 x 10° 
Senne magnetized 0.00143 0.00211 0.00345 1.92x10° 2.83x10° 4.63 x 10° 
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Fig.8 The average acceleration a, the relative escape rate 3, and the effective collision frequency v/;, of three chaos regions 


versus the chaos region radius R,/Rp in the double Y-type magnetic configuration! 


在 模拟 中 , 选择 2 x 105 个 粒子 , 3 


这 些 粒子 初始 


均匀 分 布 在 -1.0 < x y < 1.0. 2’ = 0 的 区 域 , 并 


始 


近 


对 
中 


具有 热 速度 为 vr = 0.3vA 的 各 向 
速度 分 布 . 此 外 , 选择 X 型 中 性 
区 半径 R。 = 2.0R,. 图 10 
HÆR Je = 0.4 和 e = 0.8 时 , 平均 加 速度 a、 相 
J aa 频率 fr 随 引 
别 对 应 了 XX 型 


的 混沌 


逃逸 率 9 和 
两 子 图 的 上 下 行 


同性 麦克 斯 韦 初 
点 和 O 型 中 性 点 附 
0 (a) 和 10 (b) 分 


导 场 的 变化 , 其 
1 区 域 和 0 型 区 


域 下 的 结果 . 从 ee 发 现 所 有 情况 下 平均 加 
速度 a 都 随 着 加 速 电场 MA 


Hr 
THY 


引导 场 存在 时 , 带电 粒子 会 在 混沌 


电 粒 子 的 加 速 主要 是 由 


的 增加 而 增 大 , 这 说 明 
加 速 电 场 决 定 的 . 当 强 
区 待 更 长 的 时 
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间 , 从 而 使 带电 粒子 的 加 速 更 有 效 . 然而 , 在 X 型 
区 域 下 对 于 弱 引 导 场 B。< Bo, 平均 加 速度 a 随 


导 场 Bs/ Bo 表现 出 复杂 的 变化 . 图 


的 相对 逃逸 率 6 随 引 


小 , 说 明 强 引导 场 会 有 效 阻碍 粒子 从 O 型 
区 域 中 当 引 导 


逃逸 . 而 在 X 型 


导 场 B,/Bo 的 增 大 而 急 
! 区域 


对 逃逸 率 8 随 引导 


场 MA = 0.001 和 0.0005. 当 引 导 场 B。< 0.5Bo 时 
相似 于 平均 加 速度 a、 相 对 逃逸 率 6 在 义 型 

引导 场 B, 表 现 出 复杂 的 变化 , 这 直接 导致 了 X 型 
区 域 的 有 效 碰撞 频率 wr 随 引导 场 表现 出 复杂 的 
况 下 的 有 效 


变化 . 在 0 型 区 


又 域 中 ， 


B, > 0.5Bo 时 


， 相 


中 O 型 区 域 
F yk 


g 


导 场 的 变化 不 大 , 尤其 对 于 加 速 电 


区 域 随 
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1J 


对 于 se = 0.4 情 
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碰撞 频率 as 随 引 导 场 的 增 大 基本 保持 不 变 , 而 对 在 X 型 区 域 中 随 e 的 增加 而 变化 不 大 , 说 明 在 磁 岛 内 
Fe = 0.8 的 情况 vg 随 引 导 场 的 增 大 而 减 小 . 通 。 ”党 引导 场 方向 的 加 速 可 以 被 强 磁 岛 阻 碍 . 从 图 11 中 ， 
过 对 比 图 10 (a) 和 图 10 (b), 发 现 除 了 引导 场 B。< 还 可 以 发 现 相 对 逃逸 率 8 和 有 效 碰撞 频率 vq 在 O 型 
0.5Bo 的 情况 下 , 有 效 碰撞 频率 war,， 即 混沌 感应 电 区 域 中 , RE Sa HE Rye < 0.5 的 情况 随 磁 岛 参数 e 增 
阻 7eg 在 X 型 区 域 中 的 值 一 般 要 大 于 O 型 区 域 的 值 . 加 而 稍微 增 大 , 对 于 强 磁 岛 s > 0.5 的 情况 随 e 的 增 


— 


在 弱 引 导 场 和 低 磁 岛 参数 es = 0.4 的 情况 下 , KEK 大 而 快速 减 小 . 此 外 , 图 11 显 示 了 有 效 碰撞 频率 za， 
域 的 有 效 碰撞 频率 vg 接近 或 小 于 O 型 区 域 的 值 . 即 混沌 感应 电阻 re 在 弱 磁 岛 s < 0.5 情 况 下 , X 型 区 
给 定 引 导 场 B。 = 1.0Bo 和 加 速 电 场 MA = 域 中 的 值 要 比 O 型 区 域 中 的 值 高 1 个 量 级 ,而 在 强 


0.001, 图 11 显 示 了 平均 加 速度 w、 相 对 逃逸 率 5 和 Bi the > 0.5 情 况 下 , X 型 区 域 中 的 值 要 比 O 型 区 域 
效 碰撞 频率 vg 随 磁 岛 参数 e 的 变化 . 发 现 平均 加 中 的 值 高 2 个 量 级 . 
速度 a 在 O 型 区 域 随 磁 岛 参数 e 增 加 而 明显 增 大 , 但 


EN a 
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> K DD 


图 9 HERSE = 0.4 时 , DEE PXE Dx Soh EMOR EX Sak FS) 


Fig.9 The X-type (top) and O-type (bottom) regions in the magnetic island-chain structure, where the magnetic island 


parameter ¢ = 0.4 has been used 38 


以 上 结果 说 明 混 沌 感应 电阻 在 无 碰撞 重 联 电 现 混沌 感应 电阻 与 经 典 碰撞 电阻 的 比值 在 X 型 区 域 
流 片 中 扮演 了 重要 角色 , 在 典型 空间 等 离子 体 环 约 是 107, 在 O 型 区 域 的 比值 约 是 104-105, 这 符合 太 
境 中 , 混沌 感应 电阻 可 能 远大 于 经 典 碰 撞 电 阻 . 例 EE DE HP sc ey Fa, BEL EK BERS) 表 2 的 最 后 一 列 给 出 
如 日 冕 环境 下 的 典型 等 离子 体 参 数 为 磁场 Bu ~ 了 X 型 区 域 下 混沌 感应 电阻 与 O 型 区 域 下 混沌 感应 
500 Gs、 等 离子 体 密度 no ~ 10° cm 习 和 温度 T ~ 电阻 的 比值 xog Noor, 可 以 发 现 X 型 区 域 下 的 混沌 
200 eV. 表 2 列 出 了 不 同 磁 岛 参数 ec< 下 X 型 区 域 和 O 感应 电阻 要 比 O 型 区 域 下 的 混沌 感应 电阻 高 1-2 个 
型 区 域 中 的 有 效 碰撞 频率 .gqg 和 vbg, 并 分 别 给 出 数量 级 , 这 说 明 混 沌 感应 电阻 在 具有 岛 链 结构 的 重 
了 X 型 区 域 和 O 型 区 域 下 的 混沌 感应 电阻 与 经 典 碰 联 电流 片 中 分 布 不 均匀 , 主要 集中 在 重 联 电流 片 
撞 电 阻 的 比值 xeg /wxcal 和 woeg /Woca. 从 表 2 中 , 发 的 X 型 区 域 . 
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出 了 X 型 区 域 混沌 感应 电阻 和 O 型 区 域 混沌 感应 电阻 的 比值 


resistivity in the X-type region to that in the O-type region 
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Fig.11 The average acceleration a, the relative escape rate 3, and the effective collision frequency Visp in the X-type and 


表 2 不 同 磁 岛 参数 下 , X 型 和 OO 型 区 域 中 的 有 效 碰撞 频率 , FA SP Se ee BBA. 表 中 最 后 一 列 给 


Table 2 The effective collision frequency of the X-type and O-type regions and the ratio of the 
chaos-induced resistivity to the classical collision resistivity in the corona are given under different 


parameters of the magnetic island. The last column in the table gives the ratio of chaos-induced 


感应 电阻 的 定量 计 

过 模拟 研究 发 现 混沌 
电流 片 中 反常 电阻 的 产 
在 2017 年 ， 将 


\ 有 非 零 引 


FA VE, a 


现 混沌 感 ) DA 电阻 i 


感应 电阻 在 


E 重 联 区 


下 的 混沌 感 
经 典 碰 撞 电 


[一 


jek 


X 中 的 分 布 是 不 均匀 


ChinaXiv 合 作 


HAFI 


64 卷 E 桢 等 : 无 磁 撞 重 联 区 中 的 混沌 感应 电阻 3 其 
的 , 这 非常 符合 反常 电阻 的 要 求 B7-381. 然而 , 目前 21] Biskamp D. PhF1, 1986, 29: 1520 
只 研究 了 国定 电磁 场 中 的 混沌 感应 电阻 , 而 空间 真 22] Lin S C, Liu Y H, Li X C. PhP1, 2021, 28: 072109 
实 等 离子 体 环境 下 的 电磁 场 要 复杂 得 多 . 如 果 磁 场 23] Galeev A A, Sudan R N. Basic Plasma Physics II. Vol- 
h i : ae hau x ume 2 of Handbook of Plasma Physics. Amsterdam: 
重 联 中 存在 扰动 场 , 会 对 单个 粒子 的 运动 轨迹 有 哪 North-Holland Physics, 1984: 271 
些 影响 ? 混沌 感应 电 阻 是 否 会 增 大 ? 如 果 使 用 全 粒 24] Bychenkov V Y, Silin V P, Uryupin S A. PhR, 1988, 
子 模拟 , 混沌 感应 电阻 能 够 为 磁场 重 联 所 需 的 反常 164: 119 
| 阻 提供 多 少 贡 THR? 这 些 都 是 尚未 解决 的 问题 ATi 由 25] Gardner J P, Katz N, Hernquist L, et al. ApJ, 2003, 
机 bea ee O 2 587: 1 
此 可 见 , 无 磁 撞 磁场 重 联 中 的 混沌 感应 电阻 值得 更 26] Petkaki P, Freeman M P. ApJ, 2008, 686: 686 
多 的 关注 和 研究 . 27] Voitenko Y M. SoPh, 1995, 161: 197 
ee 28] Bellan P M. PhRvL, 1999, 83: 4768 
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Chaos-induced Resistivity in Collisionless Magnetic 
Reconnection Region 


WANG Zhenl23 CHEN Ling? WU De-jin? 


(1 Institute of Space Physics, Luoyang Normal University, Luoyang 471934) 
(2 Purple Mountain Observatory, Chinese Academy of Sciences, Nanjing 210023) 
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ABsTRAcT Collisionless magnetic reconnection, which converts the magnetic energy into the kinetic 
energy of plasma particles via the heating or acceleration, has been believed widely to be able to explain 
various eruptive phenomena such as solar flares and geomagnetic storms. However, the microphysical 
mechanism of anomalous resistivity in the collisionless magnetic reconnection is still an unsolved funda- 
mental problem. Among the many physical mechanisms of anomalous resistivity formation, chaos-induced 
resistivity based on the chaos of the charged particle orbits near the magnetic neutral point is not the most 
popular formation mechanism, but its microscopic physical picture is the clearest. This paper first briefly 
reviews the early research and physical model of the chaos-induced resistivity in collisionless magnetic re- 
connection region, introduces the recent research progress of the chaos-induced resistivity, and expounds 
the future research direction of the chaos-induced resistivity. 


Key words magnetic reconnection, anomalous resistivity, chaos-induced resistivity, Sun: coronal mass 
ejections (CMEs) 
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